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Abstract 

We discuss the front propagation in the A + B — > 2A reaction under subdiffusion which is 
described by continuous time random walks with a heavy-tailed power law waiting time probability 
density function. Using a crossover argument, we discuss the two scaling regimes of the front 
propagation: an intermediate asymptotic regime given by the front solution of the corresponding 
continuous equation, and the final asymptotics, which is fluctuation-dominated and therefore lays 
out of reach of the continuous scheme. We moreover show that the continuous reaction subdiffusion 
equation indeed possesses a front solution that decelerates and becomes narrow in the course of 
time. This continuous description breaks down for larger times when the front gets atomically 
sharp. We show that the velocity of such fronts decays in time faster than in the continuous 
regime. 

PACS numbers: 
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I. INTRODUCTION 



Reactions under subdiffusion have attracted much attention in recent years due to their 
growing relevance for description of processes taking place in porous media such as certain 
geological formations or gels, in the crowded cell interiors and in many other strongly in- 
homogeneous environments including modern drug delivery systems. We focus here on the 
autocatalytic conversion A + B — )■ 2A, a reaction that exhibits travelling front solutions 
if the initial conditions are chosen appropriately, i.e. if A and B are initially separated in 
space jl, 2|. 

We concentrate on situations when subdiffusion can be modelled within the CTRW 
scheme with a waiting time probability density function (pdf ) decaying according to a power 
law, ij}(t) oc t^ x ~ a . The continuous description of the A + B — > 2A reaction under subdiffu- 
sion. following locally the mass action law corresponding to the FKPP equation, was derived 
in [3( and is given by a partial integro-differential equation with a kernel depending on the 
particle concentrations at all times. 

In that preceding work we have shown analytically that the resultant minimal front 
velocity goes to zero under the assumption of a constant front shape, which was interpreted as 
propagation failure. In a following paper j^jj], numerical simulations corroborated this picture, 
while two different regimes of front propagation were identified. In the fluctuation dominated 
regime, pertinent to large reaction rates, the front velocity was found to decay as v(t) oc t a ~ l , 
whereas in the regime of small reaction rates, for which the continuous description applies, 

a — 1 

the front velocity was observed to go as v(t) oc . Longer simulation runs of the continuous 
case (small reaction rates) revealed that after an intermediate regime that ranged over less 
than two orders of magnitude in time where v(t) oc t 3 ? applies, the exponent sets in to 
decay [5(. Hence the alleged exponent conjectured from the continuous picture was not the 
final one. Up to now, there has not been any physically sound interpretation of the front 
velocities found in these simulations. 

In this work we attempt to fill this gap by giving a crossover argument that is used to 
construct an Ansatz for the solution of the reaction subdiffusion equation at the leading 

a — 1 

edge. We found that in order to maintain a front velocity that goes as v(t) oc , the 

a.— 1 

additional assumption of the width of the front going as has to be made, so that 
the front does not maintain a constant form in the course of its propagation. Since the 
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front's width decreases with time, any real (or simulated) subdiffusive FKPP system will 
sooner or later undergo a change of regime: the front will get atomically narrow and the 
continuous scheme breaks down. Physically this has to do with the fact that at large times 
the jump rate always becomes small compared to the reaction rate, so that the fluctuation 
dominated regime sets in. Since the particles react before they are able to leave the site, 
the front becomes atomically sharp. We suggest that the findings in (decay of the 
exponent characterizing the time dependence of the front velocity) can be interpreted in 
the sense of a transition from the intermediate asymptotics of the reaction described by 
the continuous reaction-subdiffusion scheme to the final asymptotics corresponding to the 
fluctuation dominated regime. We start by presenting simple physical arguments in favor 
of this picture. We then show that the intermediate asymptotics with v(t) oc t~^r indeed 
appears as a possible solution of the corresponding integro-differential reaction-subdiffusion 
equation. Physical arguments show however that this asymptotics cannot be the final one, 
and that the final regime is fluctuation-dominated. We then turn to a numerical investigation 
of this fluctuation-dominated regime and show that the subdiffusive nature of the motion 
leads to additional fluctuation effects absent in the normal diffusive case. 



II. CROSSOVER ARGUMENTS 



Under normal diffusion and with the overall particle concentration A+B = c being locally 
conserved, the A + B — > 2A reaction is described by the Fisher-Kolmogorov-Petrovskii- 
Piscounov (FKPP) equation 

dA ^ = DAA(x, t) + k(c - A) A 

that has been extensively studied in the past. According to its classical solution , fronts 
propagating with velocities v > ylkcD are possible, and it is moreover known that for step- 



like initial condition the solution with minimum speed, v = V^kcD, is the one which is 
really achieved at long times. 

In order to gain intuition about the front behavior under subdiffusion, we make use of 
the following idea: for any waiting time pdf ip with finite mean (t), the behavior at very long 
times t 3> (t) corresponds to normal diffusion, so that the behavior pertinent to reaction- 
diffusion schemes is recovered only if time t is large enough. On the other hand, if the 



3 



initial domain of the pdf can be approximated by a power-law, ip(t) oc t -1- " up to some 
truncation time T, the behavior at short times should correspond to the one in subdiffusion, 
and there must be a smooth crossover from one regime to the other. We therefore consider 
the truncated power-law waiting time distribution with truncation parameter T, 

,a m - ito + Tr at ° err t) m 

with mean value 

m aTtg + t (tg - (T + t )") 

W (a - 1) (t? - (T + to)-) ■ 1 j 

For T > t , (t) « Y^tgT 1 "". 

For small times t « J 1 , when the system does not feel the cutoff, the behavior of the 
velocity will be similar to that in subdiffusion, whereas for large times the behavior will be 
the classical one with a constant minimal velocity. The crossover between the two regimes 
must thus take place at some crossover time t cr . We assume that in the anomalous domain 
v oc t' 3 , and that after this a crossover to normal behavior sets in. In the case when the 
normal behavior is described by the FKPP scheme this corresponds to v — const. ~ V ckD, 
with D = a 2 /2(t), where a is the step's length of the corresponding random walk process (an 
irrelevant microscopical variable), and the time behavior of the velocity in the anomalous 
regime is given by the equation 

.2 1 1/2 

(3) 



L cr — 



ck- 



2(t(t cr ))_ 

In order to determine the crossover time we concentrate on the most basic quantity that 
is known in the normal as well as in the anomalous case, i.e. the number of performed steps, 
a measure of mobility, which is given by 



n D (t) = ^ (4) 



in the normal regime t t cr , and 



t a 

n SD(t) = M 7- (5) 
I [1 + a\t'o 



in the subdiffusive regime t t cr . 

By enforcing nsD(t cr ) = n£)(t cr ) we find 



and hence t cr oc T (more precisely t\~ a = r( . 1+ct ^ 1 _ a ^ T 1 ~ Q! ). Obviously, the larger we choose 
the cutoff-parameter T, the larger becomes the crossover time. At the time the crossover 
takes place, the quantities characterizing the behavior of the system, such as the number of 
performed steps, the front velocities etc. have to match for the two regimes. Tuning T we 
get the respective values of the quantities of interest at t cr , for example the mean waiting 
time (t) oc t x ~ a for the normal case in terms of t cr . From Eq.(j3]) we then get 

v(t < t cr ) oc t 3 ^ 1 (7) 



in the subdiffusive regime. Correspondingly we can define other time-dependent effective 
characteristics in the anomalous regime, e.g. an effective mean waiting time, (t) e ff oc 
(the parameter t cr is changed to t) which yields an effective, time dependent diffusion 
coefficient -D e // oc from which Eq. (f7|) can be obtained via the classical formula 

v = ^2ckD e ff. This discussion elucidates the source of the anomalous front velocity in the 
regime of small reaction rates, as found numerically in J4]. 

We note that even the case for normal diffusion is not simple at all, especially when 
the one- dimensional situation is considered, the one especially prone to fluctuation effects. 
To understand the situation we first recall that the FKPP equation, if it holds, has the 
same form in whatever spatial dimension, and provides us not only with the velocity of 
the front, but also with the front's width. Since in any spatial dimension d the dimensions 
of the concentration [c] = L~ d and that of the reaction rate [k] = T _1 [c] _1 are connected 
to each other, so that [kc] always has the dimension of the inverse time, the combination 
V Dkc always has the dimension of velocity, and the combination w = \J D/kc always has 
the dimension of length. The characteristic width of the front is thus proportional to our 
parameter w, see Q for a quantitative discussion. The velocity of the front and its width w 
are connected by a simple relation 

v ~ w/r = D/w (8) 

where r = w 2 / D is of the order of the time which it takes a particle to diffuse through the 
front's width. 

Here it is important to note, that the width w is the only relevant parameter of the 
dimension of length in the continuous theory, but going to the particle picture, another 
characteristic length, the interparticle distance I = c~ l l d emerges, and an additional dimen- 
sionless parameter U = w/l appears. The parameter II gives us the front width measured 



in the units of the interparticle distance, and quantifies the strengths of fluctuation effects 
in the A + B — > 2B reaction. 

According to the Buckingham's H-theorem, the velocity v, [v] = LT _1 has to depend on 
the parameters of the problem as 



with / being a yet unknown function of a dimensionless parameter, and the prefactor of 
/ reproducing the classical FKPP behavior of the velocity. The prefactor of / has the 
same form in any spatial dimension, while the dimensionless argument of the function / has 
different form in spaces of different dimension. Evidently, the continuous description only 



works if w ^> I, i.e., in the classical case where ■sJT) Jck ^> c~ l l d : for large concentrations 
and diffusion coefficients and for small reaction rates. In this case there are many particles 
within the front region, and the continuous description does hold. For II ~ 1 corresponding 
to the atomically sharp front, the number of particles across the region fluctuates strongly, 
and therefore front propagation is fluctuation dominated. 



Ill and fl. 



The 



Let us now concentrate on the one-dimensional CclSC, clS discussed in 
fluctuation doufluated r e gim e iu Id correspouds to u cc Dc Q, wh.cu cau be easfly u„de r - 
stood within Eq.flSJ) by assuming the width of atomically sharp front to correspond to the 
interparticle distance, w ~ I — c _1 . Repeating the same crossover arguments, as in the 
previous case, this kind of behavior under normal diffusion is mirrored onto the form 

v(t) oc t"- 1 (9) 

for the velocity time dependence in the subdiffusive case. 

The same crossover arguments as applied to the velocity, can be also extended to the 

1/2 

width of the front. Since the front width w oc D e ^ is a decaying function of time in the 
subdiffusive case, the condition for continuous description to hold breaks down for times 
long enough, and the transition from the intermediate "classical" asymptotics, Eq.flT]), to 
the final fluctuation-dominated asymptotics Eq.([H]) inevitably takes place. 

In what follows we first show that the "classical" asymptotics, Eq.(JJJ), indeed appears as 
a possible solution of the reaction subdiffusion equation, and then we change to investigating 
the far asymptotic regime, when the reaction-subdiffusion equation breaks down. This is 
done by use of extensive numerical simulations. 



III. CONTINUOUS REACTION-SUBDIFFUSION REGIME 



Let us assume the front to behave in accordance with our crossover arguments, namely 
to have the velocity and the width going as t 2 ^ (i.e. with position x(t) oc v^t 1 ^). The 
overall form of the front will be assumed exponential at its leading edge x — > oo. Thus, the 
following Ansatz is made: 



A(x, t) = A exp 



l-a ( 1 + a 

-X t 2 [X — V t 2 



A exp 



-\ t 2 Z 



(10) 



where z = x — vot^ is the comoving variable. (The exponential Ansatz is due to the 
fact that we will anyhow linearize the equations at the front's far edge, and we know from 
elsewhere jg] that the (stationary) solutions of linear react ion- sub diffusion equations are 
exponentials.) 

The equation for the concentration of A-particles A(x, t), with c being the overall particle 
concentration, is (cf. |3j) 

dA(x,t) ~ 2 rt 



dt 



a 2 r 

k(c — A(x, t))A(x, t) H A / M(t-t') 

2 Jo 



x (A(x, t') — c) exp 



kA(x,t")dt'' 



dt'. 



(11) 



exp 



We note that A(x,t) becomes small at the leading edge x — > oo, and 

rt 



1, so that 



f t ,kA(x,t")dt" 

1^ A^M{t-t'){A{x,t') -c)exp -JfejT A{x,t")dt" | dt' 



dA(x,t) a 2 
dt ~ ~2 

+k(c- A(x,t))A(x,t) 
„i rt 



(12) 



/ M{t-t') AA{x,t') -2VA{x,t') [ kVA{x,t")dt" 
Jo Jt 1 

+(c - A(x, t')) J kAA(x, t") dt" - (c - A(x, t')) ( J WA(x, t") dt 1 



+k(c- A(x,t))A(x,t) 



dt' 
(13) 



In particular, with Ansatz ( JTUl) and taking into account that the term t ~ is negligible for 
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large t, we have 

dA(x,t) 
dt 



A exp 



— Xot 2 ( x — v t 



l + ct 
2 



1-g a-1 OL + 1 1+q 1 — Oi ( l + a 

w A t 2 £ 2 — X t 2 — - — I x — v t 2 



A exp 
v X - 



—Xot 2 ( x — t> t 
1 — cc 



l+a 



~ A exp 
Vi(x, t) = -A Xot^ exp 
AA(i,t) = A Aot 1_Q exp 



— A t 2 (a; — v Q t 



ct+l , 
2 



wo A c 



. l + a 

-A t 2 x — t> t 2 



-X t 2 x — v t 2 



Proceeding as in [3( we have to first order in concentration for the A-particles: 

dt' 



™\ j ^-^AAoexp -X t'^r(x-v t'^) 



a 



2 /■* 



+ - 



M(t - t')ck / AA exp -A t"^ (a; - v t"~ 



+ckA(x,t), 



(14) 



(15) 



i.e. 



q+l . 



A exp — A t 2 (x — w t 2 ) 



v X 



— I M(t -t')A X 2 t n - a exp -X t'—{x-v t'— 



'0 

,2 rt 



dt' 



/ M(t-t')ckA X 2 / t'^exp [-X t"^ (x- v t" 



, l + a 

2 



q+l , 



+c/cy4 exp — Xot 2 (x — w t 2 ) 



with the kernel 



M(u) 



dt" dt' 



(16) 



1 - lj){u) 

in Laplace domain (which corresponds to the Riemann-Liouville fractional derivative of order 
1 — a in the subdiffusive case, f^y^ J * ( t _ t / 1 )i- a (-)dt'). 

We note that in the following we assume ip(t) oc r a t~ 1 ~ OL so that the new parameter r 



and the old one to from the original waiting time distribution ijj(t) 



at: 



(t+t ) 



1^ (i.e. the ip we 



truncated for the crossover argumentation in the preceding section cp. (JT])) turn out to be 
the same, r — to. 
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Altogether we have then for z = x — Vot^~ and t large: 



X v Q exp 



J. — Lt 

-A t 2 z 



exp 



-A t 2 2; 



2lYa)r(l - o)t c 



B\l + B] 



+ ck 



(17) 



where B is a constant that originates from the estimation of the involved integrals, see 
Appendix [Aj with B(a,2 — a) > B > and B(i>,fi) being the Beta-function. This yields 
the dispersion relation for Aq: 



= A* + ^ 



[l-B]-Vo ck 
— A + 



K*B 



K*B 



(18) 



with 



OT , / — t — r - = = t^t, where K a is the generalized diffusion constant. From 

21 (1— a)T a a 1 (o) ' " 



ckKt 



we find the restriction 




2K*B 



AKfB 2 



K*B 



,ckK* a 
^ ^0 



[1 - B] - v ) 2 > 4ckK* a B, 



(19) 



(20) 



a quartic equation in v which yields 



1 + B±2VB 



(21) 



Note that in the normal case B — 1, the minimal front velocity v m i n = ±2\/ cDk is repro- 
duced; the other solution is a double one at v = for which there is no front. Recall again 
that B(a, 2 — a) > B > 0, therefore eq. ( 121ft always has real roots {B(a, 2 — a) > 1 for all 
a < 1). 

This analysis shows that there exists a set of (nonzero) parameters Ao and Vq for which 
Ansatz (|T0|) yields a solution to the linearized reaction subdiffusion equation ( fT3|) . although 
the integrals appearing in the calculations can only be estimated approximately. We note 
that neither an Ansatz taking a front velocity going as v (t) oc t a ~ l nor an Ansatz with v(t) oc 

a — 1 

t~^~ and a constant front width yield an asymptotic solution of the reaction-sub diffusion 
equation, and therefore such types of behavior are impossible within the continuous scheme. 

In our previous simulations we were not able to detect the changes in the front shape, 
presumably due to our averaging procedure over several runs, and hence did not conjecture 



any change of regime in 



4]. On the other hand, our simulations were not carried out for long 



enough times to detect the change of regime in the velocity variable. Since such transitions 



take place only very slowly, muc 



more extensive simulations of [5] showed. This suggests that indeed the continuous regime 
as considered above does not describe the final behavior of the front. Now we can interpret 

n 

the findings of Ref. |5j as the setting in of a slow transition to the fluctuation dominated 
regime. 



i longer runs of the simulation were really necessary, as the 



IV. FAILURE OF THE CONTINUOUS DESCRIPTION: ATOMIC ALLY SHARP 
FRONTS IN SIMULATIONS FOR LARGE TIMES 

Since the subdiffusive front is slowing down and becoming steeper in the course of time, 
any system will sooner or later enter a regime already discussed in ^| for sub diffusion and 



y syst 

flfl 



in p, [7| for normal diffusion. This regime is a fluctuation dominated one and is no longer 
described by continuous approaches. Since the integral kernel M(t) of the linearized reaction- 
subdiffusion equations decays and determines the mean density of steps in time, the waiting 
times for particles at a site become so large in the course of time that the motion of the 
front is governed by the first A-particle entering a new site. All B-particles at the same 
site have enough time to react with A before the next jump from the site takes place, the 
reaction rate dependence disappears, and the behavior of the front gets to be the same as 
in the reaction on the first contact. 

Under such a condition the velocity of the front's motion can be estimated using the fol- 
lowing a^ent (adapted from flQ for onr serial updating sehe-ne). Let Us eon S1 de r 
the front position as fixed by the rightmost A-particle(s), and concentrate on the next jump 
of the front particle. If the A-particle is alone at its front position, this next jump takes place 
with probability 1/2 by an amount ±a, so that the net front displacement after such a step 
is zero on average. On the contrary, if there is more than one particle at the front position 
(the probability of which is ac if the concentration is defined as a number of particles per 
unit length) the front moves by a to the right, if the particle makes a step forward (which 
happens with probability 1/2), and does not move, if it jumps backwards, since then there 
is at least one other particle, which keeps the front position where it was. Therefore, at a 
step of a front particle, the front moves on average by a distance a 2 c/2. Since the rate at 
which the particle moves is defined by the time-integral of the memory kernel M, the front's 
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velocity is given by 

v « ^ / M(t-t')dt'. (22) 

Let us first derive the asymptotic jump rate of the particles. Consider the generic waiting 
time pdfs with the asymptotic behavior 

i/>(t) oc r a r l - a . (23) 

The (cumulative) probability to make a step until t, for t large is then 

#(t) ~ 1 - r a r a ; (24) 

or in Laplace domain, using the Tauberian theorem 

*(u) ~ --r(l-«)r a M" 1+a , (25) 

■u 

so that the pdf 

^(u) ~ 1 -T(l -a)T a u°. (26) 
The rate for a particle to jump is £ M(t - t')df or in Laplace domain 

mm = m ^ i u - a (27) 

u l-^(u) r«r(l-a) V ; 

for u — y so that we have an expression for the velocity in the Laplace domain given by 

C { v (t)\ = — t~ -u' a . (28) 

1 WJ 2 r a r(l -a) K ' 

Transforming back to the time domain yields 

= wi^ftr 1 = K °wr = cA » r "- (29) 

With r(a)r(i-a) = Sm !r 7r ' > ^ ne f ron t velocity is better expressed as 

T" Z 7T 

which corresponds to the position of the front going as 

/ x N A f f , . , a 2 sin (a7r) „ 
*(t) = ^ = l «<*)* -—-L-L*. (31) 

(A^4 is the total amount of A-particles). 
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Note that the definition of the characteristic waiting time r adopted here does not allow 
for simply taking a = 1 to perform the limiting transition to normal diffusion, as found e.g. 
for the exponential distribution of waiting times, if)(t) = (t)^ 1 exp(— tj (£)). This is due to the 
presence of the divergent T(l — a) in Eq. f[2"6"j) : the case a = 1 corresponds, strictly speaking 
to still (logarithmically) divergent mean waiting times. For the normal case with converging 
mean Eq. (1261) reads if>(u) — 1 — (t)u, and, after performing the same steps as above, the front 
velocity of the normal fluctuation dominated regime, Vfi uc t = cD, with D being the diffusion 
constant, is recovered js]. Fig. [TJ shows the total number of particles in the simulation for 
the fluctuation dominated regime, i.e. reaction on contact, for a concentration c = 0.3. In 
these simulations we had 7 runs for a = 0.9, 18 for a = 0.8, 41 for a = 0.75, 13 for a = 0.7 
and 18 for a = 0.6. 




lxlO 4 2xl0 4 5xl0 4 lxlO 5 2xl0 5 5xl0 5 lxlO 6 2xl0 6 
t 



FIG. 1: Front position for a = 0.9, 0.8, 0.75, 0.7, 0.6 (upper to lower graphs), c = 0.3. 
Red lines denote fits of the large time behavior. 



Table [J shows the exponents of the long time fits A0i = Ft 13 which coincide well with a. 



a 


0.6 


0.7 


0.75 


0.8 


0.9 




0.603 ± 0.004 


0.708 ± 0.004 


0.750 ±0.001 


0.775 ± 0.002 


0.890 ± 0.009 



TABLE I: Exponents for the fit Na = Ft 13 for different a. 



The values of the prefactor F found from the simulations turned out to be however 
larger than the predicted ones in (IHTj) by around 30 — 40%. In order to find out about 
the origin of this difference, we performed simultaneous simulations of subdiffusion and of 
subdiffusion with randomized particles, i.e. in the situation when the particles lost their 
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individual memory and were chosen randomly to jump when a jumping time was reached. 
This variant of the reaction closely mimics the behavior assumed to derive Eq. (1221) . namely 
the assumption that the rate at which the steps of the rightmost A particle are made is 
equal to the mean jump rate of all particles at time t: we fully disregard the fact that the 
rightmost A is a very special particle, with its special prehistory. 



a=0.75 




1 2 5 10 20 50 

10 4 xt 



FIG. 2: Time dependence of the total amount of A-particles Na for the subdiffusive case 
(squares) and subdiffusion with randomized particles (circles). The black line denotes the 

theoretical curve according to ( 13TT) . The inset shows the situation for an exponential 
waiting time pdf (with mean 1), t goes from 10 to 5 x 10 4 , Na goes from 6 to 2000. The 
black line denotes again the theory, Na = Dc 2 t; c = 0.3. 



Fig. [2] shows the time dependence of the overall amount of A-particles for a = 0.75. 
The theoretical curve (|3T]) lies much closer to the simulation results of subdiffusion with 
randomized particles. The remaining difference between the simulation of the randomized 
particles and the theoretical result is presumably due to the fact that convergence to the 
asymptotic behavior in subdiffusion is very slow. Apparently, the full subdiffusive picture 
implies an additional fluctuation effect. For a better interpretation of the results, we also 
simulated the case of normal diffusion. The inset of the figure shows the situation for an 
exponential waiting time pdf with mean 1, where the simulated front behavior converges to 
the predicted behavior indicated by the black line, Na = Dc 2 t. We note that Warren et 
al. |9] detected a fluctuation effect in the normal case that occurs at small concentrations. 
However, as the inset shows, due to the sequential update in our simulations, this effect 
does not come into play here and our theoretical approach is sufficient to explain the front 
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behavior in the normal case. 



a=0.75 




1.0 1.5 2.0 3.0 5.0 7.0 10.0 15.0 

10 6 xn 



FIG. 3: Total amount of A-particles for the normal case (triangles), the subdiffusive 
case (squares) and subdiffusion with randomized jumps (circles), both a = 0.75, depending 

on the total number of performed steps; c = 0.3. 

Fig. [3] shows the dependence of the overall amount of A-particles on the total amount 
of steps performed for a = 0.75. Comparing the two subdiffusive prescriptions (original 
and randomized) as well as the normal diffusion reveals that the randomized version of 
subdiffusive front behavior is more akin to the normal diffusive front behavior than the full 
subdiffusive version: If we interpret the number of steps n as the internal, operational time of 
the process, the randomized subdiffusive setting and the normal diffusive one have the same 
asymptotics, whereas the full original subdiffusive front position differs by a certain factor. 
Fig. H] shows the quotient of the original subdiffusive front position and the randomized one, 
which can be used to quantify this effect that turns out to be around at least 20 — 30%. 

Obviously, the additional fluctuation effect of the front behavior is genuinely due to 
subdiffusion. This effect cannot be explained within the mean-field description of the front 
behavior, but comes into play through the interaction of the particles at the front: The rate 
at which a front particle performs a jump is higher than the average jump rate of a single 
particle in the system. If the particle at the edge of the front is subject to a very long waiting 
time (which happens not often, but occasionally), other particles will outpace that particle 
and take the lead. Hence, the impact of very long waiting times in single particle dynamics 
on the front motion is considerably reduced. 
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t 

FIG. 4: Quotient of total amount of A-particles for the subdiffusive case and subdiffusion 

with randomized particles A ' SD function of time. 

r Na,rcp 

V. CONCLUSIONS 

We discussed the front motion in the A + B — > 2A reaction under subdiffusion described 
by continuous time random walks where the reaction is governed by the mass action law 
on a microscopic scale. We have shown that at intermediate times, as long as the process 

a — 1 

can be described within a continuous picture, the front velocity goes as v (t) oc t~5~. The 
decay of the front velocity goes along with a decay of the width of the front, which at longer 
times therefore gets atomically sharp. At such times the continuous picture, implied by 
the description within the reaction-sub diffusion equations scheme, inevitably breaks down. 
The typical time scale of diffusion becomes very large compared to the typical time scale of 
reaction, and a crossover to the fluctuation dominated regime takes place where the front 
velocity decays faster, v (t) oc t a ~ l . This fluctuation dominated regime is the same as in the 
reaction on the first contact, and is characterized by additional fluctuation effects compared 
to the case of normal diffusion. 
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Appendix A: Evaluation of Integrals 



We investigate the integrals in expression (jT6j) term by term, from left to right and take 
into account that the constant A Q cancels. 



h 



M(t - t')\ 2 t n - a exp -\ t'^r {x - v Q t'*) 



dt' 
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1 d rt 
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— X t 2 [x — v t 



x 



r(a) dt J (t — t 



——t n - a exp 



f ± '_ / , 1+Qi \ 1 — a / l-\-a 

—X t 2 a; — v t 2 + A i 2 x - t> t 2 



(Al) 



This expression can be estimated from above since t' < t: 
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-A t 2 [x — V t 2 
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1 d 



1 



T(a)dt J (l-t') 1 - 
' B(a,2-a) 



t n ~ a dt' 



-t n - a dt' 



r(a) 



(A2) 



the integral in ( 1A1I) is monotonic, i.e. it must tend to a constant value B < B (a, 2 — a) for 
large times (B = 1 for the normal diffusive case, in particular). 
We used here the definition of the Beta-function B(fx, v) = ^[^^ ■ 



J 2 = f M(t - t')ck\l [ t" l ~ a exp 
Jo Jt' 

At the far edge of the front, our comoving variable z = x — vot~^~ is very large. The 
transition to large z can be achieved by introducing a large parameter 7, so that the integral 
appearing in the integrand of (1A3|) obtains the form of a Laplace integral which allows for 
an asymptotic estimation for 7 — > 00: 
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(A4) 
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that means that for large 7 the value of the above integral is asymptotically determined by 



the points where the exponent in the integrand attains its maximum, see e.g. 
Hence, (IA3D becomes 



111. 
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v r(l — q)t 
with B < B(a,2-a), cf. (TA2 
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